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We consider cosmological constraints arising from the background expansion history on the ef- 
fective field theory of cosmic acceleration, a theoretical framework that allows for a unified way to 
classify both models of dark energy and modified gravity within the linear regime. In the Einstein 
frame, the most general action for the background can be written in terms of a canonical scalar field 
which is non-minimally coupled to matter. The leading corrections to the action are expressible 
through a quartic kinetic term, and scalar couplings to a Gauss-Bonnet curvature term and the 
Einstein tensor. We determine the implications of the terms in this general action for the predicted 
expansion history in the context of dynamical attractors. We find that each modifies the matter 
dominated and/or accelerative eras in ways that allow us to place cosmological constraints on them. 
We present current constraints on the effective action using the latest Type la supernovae, Cosmic 
Microwave Background, and Baryonic Acoustic Oscillation data. This includes finding that the 
scalar field EFT with a coupled Gauss-Bonnet term and the data are significantly discrepant. 



I. INTRODUCTION 

Current measurements of the expansion history have 
provided complementary evidence for the universe's re- 
cent transition to accelerative expansion. The primary 
constraints have come from geometric measurements of 
cosmic distances from Type la supernovae luminosity dis- 
tances [1-7], the angular scale of the sound horizon at 
last scattering imprinted in the Cosmic Microwave Back- 
ground (CMB) temperature correlations [8-13], and in 
the distribution of large scale structure through Baryon 
Acoustic Oscillations (BAO) [14-16] 

Understanding the underlying cause of this acceler- 
ation, under the broad banner of "dark energy", rep- 
resents one of the primary challenges facing contempo- 
rary astronomy and astrophysics. Though current con- 
straints are entirely consistent with a cosmological con- 
stant, A, a dynamically evolving dark energy is by no 
means ruled out by the data. Theoretically, the coinci- 
dence and fine-tuning problems associated with A have 
led theorists to consider broader alternative causes for 
the acceleration. These include introducing a new type of 
matter, typically in the form of a scalar field, or modifica- 
tions to General Relativity on large scales. Many models 
have been suggested in the literature, e.g. quintessence, 
non-minimal couplings between quintessence and matter, 
k-essence, /(i?)-gravity, Brans-Dicke theories, Galileon 
gravity, Dvali-Gabadadze-Porrati (DGP) gravity, see, for 
example, [17] for a review of a wide range of dark energy 
models, and [18, 19] for recent discussions on modified 
gravity models. 

There has been significant progress in characterizing 
possible theoretical mechanisms to generate cosmic ac- 
celeration. Arguably, however, the improvements in the 
precision and breadth of observational constraints, of the 
expansion history and the growth of large scale structure, 
are providing remarkable opportunities to learn about 
the phenomenological properties of dark energy. As the 
accuracy of cosmological measurements continues to im- 



prove at a rapid pace, it is crucial to establish a robust 
and systematic way to connect data with fundamental 
theory. In particle physics and condensed matter sys- 
tems the framework of Effective Field Theory (EFT) has 
proven to be very successful at this endeavor (for reviews 
sec [20, 21]). The procedure is to construct the most 
general theory for observables that is compatible with 
the expected symmetries of the theory. The terms in 
the action can then be constrained by a combination of 
data from experiments, as well as ensuring theoretical 
self consistency. 

The application of the EFT approach to cosmology 
has recently been carried out for single and multi-field 
inflation [22, 23]. In these papers the authors construct 
the most general action for the cosmological fluctuations 
around an assumed inflationary background. For infla- 
tion, this is reasonable, since apart from establishing a 
successful period of inflation the key observables are the 
correlation functions of scalar and metric perturbations. 
However, the same is not true for present day accelera- 
tion where the background expansion itself gives rise to 
observables, motivating the construction of the EFT of 
the background. Such an approach was first developed 1 
for inflation by Weinberg [25] and for late-time cosmic 
acceleration in [26] - with later work appearing in [27- 
30]. 

The goal of these approaches is to systematically cat- 
egorize all proposals for cosmic acceleration, other than 
a cosmological constant, so they can be efficiently scru- 
tinized by theory and experiment simultaneously. To fa- 
cilitate such an approach a Lagrangian was constructed 
that could not only account for dark energy models, such 
as quintessence, but also theories of modified gravity in 



For a review of earlier work utilizing EFT methods to address 
the importance of high energy signatures on inflation we refer 
the reader to [24] and references within. 
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the linear regime. The relevant Lagrangian is that of 
a scalar-tensor theory non-minimally coupled to gravity, 
where in the modified gravity case the scalar can be in- 
terpreted as the extra degree of freedom arising from the 
longitudinal component of the graviton. 

In [31], the authors restricted their attention to 
quintessence and utilized the EFT of the perturbations 
around the accelerating background, as was done for in- 
flation in [22]. This has the advantage that perturba- 
tions about non-linear backgrounds can be considered, 
but the disadvantage that the acceleration of the back- 
ground must be assumed a priori. Regardless, this type 
of approach is important when considering models be- 
yond the linear regime, and in particular for models of 
modified gravity where recovering the predictions of Gen- 
eral Relativity and consistency with solar system tests 
requires the EFT of the background to ultimately fail 
- e.g. in regions of high density leading to a screening 
mechanism e.g. [32-34]. 

In addition to the EFT background approach, one may 
also consider imposing symmetries on non-linear back- 
grounds to retain a well defined Cauchy problem. In the 
EFT approach, this is never a concern since for energies 
below the cutoff of the theory, the number of degrees 
of freedom remains fixed, and higher time derivatives 
never appear in the equations of motion [25] . However, as 
shown long ago by Horndcski [35] , by restricting the op- 
erators to be considered in the action it is possible to con- 
struct non-linear backgrounds that give rise to only two 
derivatives acting on fields at the level of the equations 
of motion. These models are less general than the EFT 
approach, but in addition to being able to capture non- 
linearities they also exhibit interesting self-tuning prop- 
erties. The authors of [36] recently studied these actions 
identifying four terms that are important for classifying 
these scalar-tensor theories. Given our EFT approach, 
these terms are captured at leading order by the action 
to be considered below, but the full non-linear effects will 
not be captured. 

A powerful approach to understanding the implications 
of a given theory is to find dynamical attractor solutions 
for cosmological evolution. These provide predictions of 
evolutionary trajectories that a theory will naturally tend 
towards, largely insensitive to assumptions about initial 
conditions. Dynamical attractor analyses have allowed 
inferences about the viability of dark energy theories, in 
the context of minimally coupled quintessence theories 
[37, 38], non-minimal couplings in the dark sector [39], 
and f(R) theories [40], and see [41] for a rather compre- 
hensive review of dynamical attractors in dark energy. 
There are ways to evade dynamical attractors, for exam- 
ple by choosing forms of self-interactions that, by con- 
struction, do not allow attractors [42] . While such selec- 
tions provide a proof of principle, they explicitly require 
a fine-tuning by the specificity of their form, which the 
generality of dynamical attractors works to avoid. 

In this paper we consider the implications for dynami- 
cal attractors in the background expansion history of an 



effective theory of cosmic acceleration expressed in terms 
of the low-energy effective action proposed by Park, Wat- 
son, and Zurek [26] and Bloomfield and Flanagan [27]. 
Our results here apply strictly to the linear regime. There 
has been recent work looking at observational implication 
of the non- linear EFT of perturbations [43-45]. 

In section II we outline the effective field theory action 
which we consider in the paper and discuss the origins 
of the contributing terms and the resulting equations of 
motion. In section III we present dynamical attractor 
solutions for the components of the effective theory and 
consider the analytical and numerical implications for the 
background expansion history. In section IV we outline 
the cosmological data sets we use to constrain the effec- 
tive theory and summarize our results. We bring together 
our findings and discuss implications for future work in 
section V. 



II. EFFECTIVE FIELD THEORY FOR DARK 
ENERGY 

Working in the Einstein frame, we consider the effec- 
tive theory for the background with leading scalar and 
gravitational corrections [26, 27] 

r { m 2 i 1 

S = J d A x^~g | -fR - -(V» 2 - V&) | 

+ J d A Xsf^, {/ g „„ rt ic(0)(V(A) 4 + fcurv^G^V^A 

+/gbW (R 2 - ABTR^ + R liVtrp R' M " rp ) } 

+s , m [n(0,V0)s„„,Vm] (i) 

where M p = = 2 A3 x 10 18 GeV is the reduced 

Planck mass and tp m are the standard model and dark 
matter fields. The first line gives the leading order terms, 
a canonical scalar field with an arbitrary potential V((f>). 
The second line gives the leading derivative corrections: 
a quartic kinetic term and a direct coupling between the 
scalar gradient and the Einstein tensor. The third line 
is a Gauss-Bonnet (GB) curvature term. The last line 
describes the non-minimally coupling between the scalar 
field and matter in the Einstein frame, with 

n(cf>,Vcj>) = e a ^ (l + f kin (<P)(V<j>) 2 ) (2) 

where a is the leading order coupling and a term pro- 
portional to fkin provides the next order correction. The 
coefficients f quart , fcurv, fcB and fkin, along with V and 
a, are all arbitrary functions of the scalar field. We note 
that this implies that the GB term is no longer a total 
derivative. 

We focus here on terms in the effective theory that 
may influence cosmic evolution in the matter dominated 
and late time accelerative eras. We do not consider 
higher order corrections to the energy-momentum ten- 
sor, T^ v = diag(— p, P, P, P) where p is the cosmic mat- 
ter density and P is the isotropic pressure that can arise 
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in the a general EFT. Specifically we neglect operators 
such as those proportional to (T^) 2 = (— p + 3P) 2 or 
T^Tpv = (p 2 + 3P 2 ), since these types of corrections 
would most likely lead to effects at all times, and be 
particularly relevant in the high density regime during 
radiation domination. 

The power of the EFT approach is that data from ob- 
servations can be used to constrain the possible values of 
the parameters, while at the same time theoretical con- 
sistency places additional constraints to eliminate regions 
of the parameter space. As an example of the latter, it 
is well known that radiative stability of the scalar field 
for quintessence puts strong constraints on the form of 
the scalar potential [46]. While the coefficients in the 
action are arbitrary functions of the scalar field, it is use- 
ful to expand these functions in powers of the cutoff of 
the effective theory. The cutoff is a dimensionful scale 
that characterizes the higher energy, microscopic physics 
that has been integrated out as our observational interest 
passes to lower energy scales. For example, for the ap- 
plication to dark energy this scale could result from inte- 
grating out electrons. We note, that given the cosmolog- 
ical background is dynamic (characterized by the Hubble 
parameter, H(t)), this can lead to several scales besides 
the masses of heavy particles going into determining the 
effective cutoff(s) of the theory - sec e.g. [47]. Lacking 
a knowledge of the high energy completion here, in this 
paper we will study the implications for each leading cor- 
rection by expanding functions with time derivatives of 
the scalar scaled by powers of the Hubble scale, whereas 
the scalar field will be suppressed by powers of the Planck 
mass. 

We emphasize that since we are considering the EFT 
of the background (not the perturbations) our expan- 
sion must always respect the hierarchical structure of the 
terms, i.e. higher order corrections must always be sub- 
leading to lower order terms. This is particularly im- 
portant to keep in mind in regards to the GB term, as 
this type of term is often considered with a large pre- 
factor, which can lead to a number of instabilities 2 . Here 
we will follow the EFT approach and treat the GB term 
as the first in a serious of higher derivative corrections. 
Indeed, the procedure followed in [26, 27] for arriving 
at the general action (1) requires that this term remain 
small for the reduction procedure that was implemented 
there to be valid. We also emphasize that the presence 
of the scalar dependent coefficient Jgb implies that the 
GB term is not purely topological and so can play an 
important role in the dynamics. 

Our cutoffs, chosen to scale with powers of the Hubble 
and Planck scales, represent the most optimistic case for 
observations, and therefore the strongest constraints for 



2 For a discussion of the pathologies associated with large GB 
terms and the difference with the EFT approach we refer the 
reader to [25]. 



theories. For larger cutoffs the observational effects of 
higher dimensional operators for dark energy and modi- 
fied gravity are known to be negligible. Models of ghost 
condensation offer an example, where there the cutoff is 
typically taken to be of order a few GeV, and so the evo- 
lution is indistinguishable from a cosmological constant 
[48] . For the EFT of quintessence, a GeV scale cutoff im- 
plies a similar story, and all models remain observation- 
ally indistinguishable from ordinary, vanilla quintessence. 
Thus, here we focus on the lowest possible cutoff con- 
sistent with the effective theory, leaving the important 
question of the UV completion to future work. 

Given these considerations we take the couplings to 
terms involving scalar gradients to be of the form: 



f quart 

fcurv 
fkin 



M 2 H 2 

H 2 

F k 
M 2 H 2 



(3) 

(4) 
(5) 



where F q ,F t 



F k are constants. Following common ex- 
amples from the literature, we assume exponential cou- 
plings and potentials for the interactions not involving 
scalar gradients, 



V = V exp -A 



e a = exp -2Q 



Jgb = F exp n 



M r 



(6) 
(7) 
(8) 



with Vo, A, Q, F and p constant. Our definitions of A, 
p and Q are chosen to be consistent with previous work, 
such as [37, 49], and Q = y/zJiC in [39]. While we con- 
sider an exponential potential in this analysis, power law 
potentials are also commonly considered, and though the 
findings differ in the details, typically they yield broadly 
comparable results in relation to the constraints on the 
coupling Q e.g. [39, 50]. 



A. The equations of motion 

We assume a Friedmann Robertson Walker (FRW) 
metric, 



ds 2 



-dt 2 + a 2 dx 2 



(9) 



where t is physical time and a is the scale factor and H = 
a/ a. Given a homogeneous scalar field, the Einstein field 
equations and energy-momentum conservation equations 
for the scalar and matter give rise to consistent general 
equations of motion given by the Friedmann equation, 

Pm{4>) + P~f + + V((j>) + Zfquart^ 



3M 2 H 2 



+9/ cur „ff V + 24<j>f GB H 3 



(10) 
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the acceleration equation, 



A. Leading order terms 



-3 M 2 H 2 




l-/< 



" ,rv M 2 



+ 2^+2 



(11) 



the scalar field equation, 



(1 + I2.f g uartj> 2 )4> + (3 + IVquarttf) H(j) + V' = 



M P 2(1- f ki n<j> 2 ) , 



-K 



quarts 



-WgbH 4 



H 
H 2 



+ 1 



~fcurv(6H 2 4> + 12HH4, + lSH 3 ^) - 3f curv H 2 <p 2 , (12) 
and the matter fluid equation 
p = -3H(p + P) 

/ ^ 3 + 2 - fen ^\p-3P), (13) 



Q 



2(1 - / fcm 2 ) 



where primes denote derivatives with respect to the scalar 
field and dots denote derivatives with respect to physical 
time, p and P are functions of 4> through Q and F^in, 
unless conformally coupled. 

While it is convenient to perform the dynamical at- 
tractor analysis in the Einstein frame, comparisons with 
observations are performed in the Jordan frame, as typi- 
cally cosmological redshifts in data assume baryonic mat- 
ter is minimally coupled, and not subject to fifth forces 
due to a direct scalar coupling. To transform into the 
Jordan frame we use the conformal transformation to re- 
move the coupling from the metric, with g^ v — £lg pu , 
dt = Q~ 1 / 2 dt, a = il _1 / 2 a and matter variables trans- 
forming as P = fl~ 2 P, and p = 0~ 2 p, where tildes de- 
note Jordan frame variables. 



III. DYNAMICAL ATTRACTOR SOLUTIONS 

In this section we use the background equations of mo- 
tion to look for the stationary 'attractor' solutions during 
the matter and late time accelerative eras. We will con- 
sider the impact of the higher order operators arising in 
the effective field theory on the cosmological dynamics 
compared with the behavior predicted by the leading or- 
der action. 



The standard dynamical solutions for a non-minimally 
coupled, canonical scalar field [37, 38, 41] are commonly 
written using dimensionless variables, see for example 
[41], 



x = 



M p H^/6' 



y = 



1 W 



Z = 



M P H ^3 



.(14) 



The Friedmann equation gives the fractional matter den- 
sity, f2 m , in terms of these variables 



ft m (a) 



Pn 



3M 2 H 2 



= l-x 2 -y 2 -z 2 . (15) 



The acceleration equation gives the effective Einstein 
frame equation of state parameter, we, 



H 3„ , 3 
-fl5 = 2< 1 + tt *> = 2 



l + x 2 -y 2 + ^z- 



,(16) 



and an effective scalar equation of state parameter, 



w E 



±z 2 

3 Z 



y 



x + y 2 



x 2 + y 



2 ' 



(17) 



The scalar and matter fluid equations give the evolu- 
tion for x, y, and z, 



dx 
din a 

dy 
din a 

dz 
dh\a 



= -y 




QQ m (a), (18) 



(19) 
(20) 



The dynamical attractors are given by the static 
solutions to these equations dx/dhia = dy/dhia = 
dz/d\na = 0. 

We can write a general expression for the effective 
Jordan equation of state, wj, during attractor-driven 
epochs, in terms of the coupling and we, using the con- 
formal transformation, 



3(1 + wj) 



3(1 + ^) - 2^6Qx 
1 - VGQx 



(21) 



where x is the fractional scalar kinetic energy component 
in the attractor regime. This relationship holds true even 
when the non-relativistic matter density, and its effect 
on the fluid equations, is negligible, and is well-defined 
as long as x ^ l/>/6Q. 

Two primary matter era attractor solutions arise: one 
dependent on the potential (MAT- A), and one wholly de- 
termined by the non- minimal coupling (MAT-Q). A sin- 
gle potential-driven late time accelerative attractor exists 
(ACC-A). 
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Attractor 


X 


y 




U>0 






RAD-A 




2\/3 
~3A~ 


4 

A^ 


1/3 


1/3 


1 (A+4Q1 
3 (A-4Q) 


RAD-null 











- 


1/3 


1/3 


MAT- A 


K 1 


x/f-Q(A-Q) 


3 Q 


Q(A-Q) 


Q 


2Q 


V 2 A -<3 


A-Q 


(A-Q) 2 A-Q 


Q(A-Q)-3 


A-Q 


A-4Q 


MAT-Q 







2Q 2 
3 


1 


2Q 2 
3 


4Q 2 
3(1-2Q 2 ) 


ACC-A 


A 

V6 




1 


A ~ 3 


- 1 - T 3 


, . A(A-2Q) 
1 3(1-QA) 



TABLE I: Table summarizing the properties of the principal dynamical attractors arising from the leading order 
terms in the effective action. They consist of two radiation era attractors, RAD-A and RAD-null, two matter era 
attractors, MAT-A and MAT-Q, and one accelerative era attractor, ACC-A. We give the values of the dimensionless 
scalar field dynamical variables, x and y, the Einsten frame fractional energy density, and equation of state, to^, 
for the scalar field, and the overall effective equation of state in the Einstein, we, and Jordan frame, w.j. 



For a minimal coupled scalar field, when Q — 0, the 
MAT-A attractor is the main cosmological scaling solu- 
tion [37, 38]. In this case, the scalar field evolves with 
the same equation of state, w, as the dominant matter 
component, with the attractive property that the radi- 
ation and matter dominated eras evolve as in ACDM. 
To satisfy this scaling solution the scalar must be sub- 
dominant, Qrf, < Q m requiring A > \/6. Nucleosynthesis 
puts a stronger lower bound on A in the RAD-A attrac- 
tor; the non-zero scalar density increases the expansion 
rate, altering the primordial abundances. For example 
[51] reported a dark energy density of < 0.09 which 
translates into the constraint A > 6.5. One cannot, how- 
ever, simultaneously realize a cosmological solution that 
flows from this matter era attractor to an accelerative 
expansion at late times, as acceleration with the ACC- 
A attractor requires A < y/2. A simple exponential dark 
energy model cannot provide a viable matter and acceler- 
ative era without an degree of freedom, such as a double 
potential or a feature in the potential [52] . 

In the non-minimally coupled scenario, an alternative 
matter era attractor, MAT-Q, is also present. The bene- 
fit of this attractor is that it can arise from the RAD-null 
attractor, evading the BBN constraints on A, and transi- 
tion to an accelerative ACC-A attractor at late times (for 
appropriate choice of A). It runs into problems observa- 
tionally, however, because during the matter era wj ^ 0. 
This was first highlighted in the context of f(R) theories 
[40], for which Q = l/\/6, w E = 1/9 and wj = 1/3 dur- 
ing the matter dominated era, showing that dynamical 
attractors for such theories are inconsistent with data. 
The MAT-A attractor also deviates from the dominant 
background equation of state when Q = 0, putting it in 
tension with data. 

The accelerative expansion, in ACC-A is enhanced as 
Q is increased for < Q < 1/A, or Q > (3 - A 2 )/A for 
A < \/2- In the case A > -\/3 the coupling Q must be 
larger than Q > 1/A. This could expand the range of A 



that can give rise to acceleration. 



B. Gauss-Bonnet term 

If uncoupled to the scalar, the Gauss-Bonnet (GB) 
term is a topological term which would play no role in 
the cosmic dynamics. If coupled, however, it will impact 
the expansion history [49, 53]. 

We introduce an extra dimensionless variable for the 
GB term 



& fGB H2 

Mr, 



(22) 



The modified Friedmann equation (10) yields a matter 
density 



fl m (a) = 1 - x 2 - y 2 - z 2 
and the acceleration equation gives 



V6. 



xv 



(23) 



2 H 
3iP 



l-v[V6: 



1 



-2,2 2 

-z + x — y 



-v Xy 2 + QVL m {a) + 2/ix 2 -v- 



(24) 



The fluid equations for y and z have the same forms as 
in (19) and (20), while the dynamical equation for the 
kinetic term, x, and the GB term are given by 



dx 
din a 



dv 
din a 




Xy + Q^m{a) -v[l 



H 

#2 



n I H 



(25) 
(26) 
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The matter dominated era allows static solutions with 
v = 0, so that the GB term is negligible and attractors 
are the same for the leading terms: MAT-A and MAT-Q 
given in Table I. The difference in including the GB term 
occurs in the late time accelerative era. An alternative 
attractor exists, which is independent of the coupling Q, 
given by 



Attractor 
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f^0 


W0 


w E 


Wj 


ACC - GB 





l 


A 


f 


-1 


-1 


-1 



(27) 



This allows cosmological evolution following the MAT-A, 
with a scaling attractor in the matter dominated era, to 
transition to a Gauss-Bonnet induced, de-Sitter solution 
at late times. Note that the accelerative era arises be- 
cause the gradient in the GB term becomes comparable, 
and opposite to that of the potential, however the energy 
density in the GB term, given by ^/6xv } tends to zero in 
both the matter and de-Sitter epochs. 

At late times v = A, y = 1 and x = 0, and <f) is con- 
stant. The time at which the transition from the matter 
dominated to accelerative era occurs depends on the rel- 
ative importance of the potential, Vq , and Gauss-Bonnet 
term, Fq, in the scalar field equation. We can obtain an 
approximate relationship between the two by assuming 
the accelerative attractor is reached today, with H = H Q , 
[53], 



A « 8/£ 



GB 



Hi 



-8F of iH 2 



Ml 



exp - 



M n 



= 1 



V 



exp 



-A- 



(28) 
(29) 



This gives an approximate relation to estimate the value 
of Fq required to give acceleration today, in terms of the 
potential, 



rpest 
t 



AM 2 (3M 2 H 2 



8^H 2 



Vn 



(30) 



This relation is only a rough estimate, as the transition 
occurs prior to today and we have not yet reached the 
pure accelerative era. It gives a sufficiently good starting- 
point, however, to guide a nuisance parameter in the 
MCMC analysis as discussed in IV A. 



The acceleration equation becomes 



H 

IP 



1 + x 2 



y 



z 2 + 12F n x 4 



with an effective, overall equation of state 

o o 1 



W E 



y 



z 2 + 12F q x 4 . 



(32) 



(33) 



The fluid equations for y and z have the same forms as 
in (19) and (20), while the dynamical equation for the 
kinetic term, x, is given by 

(1 + 12F q x 2 )x' = -xU + {l + Wx 2 F q )^+72F q x 2 \ 



-Qn m (a). 



(34) 



The RAD-null solution during the radiation dominated 
era is unchanged. For the RAD — A scaling solution z is 

changed to z = V A 4 ^ 2 56F g however this just repre- 
sents an adjustment of the relative contributions of 7 
and that keeps the effective equation of state unal- 
tered, we = 3- The same is true for the MAT—X scaling, 
the solution for y is altered but w E = j^q is unchanged. 

The coupling in the matter dominated era still admits 
the MAT-Q attractor, but with a coupling that is depen- 
dent on both Q and F q . 

Considering all the terms in x' = one finds a matter 
attractor solution, xmat, that differs from the leading 
case 



= 1296F 2 x 7 MAT 
+3(1 - 12F q )x 



U4F q x 5 MAT 



MAT ~ V®Q X MAT 



2>&V&QF q x A MAT 

3xmat + V&Q- (35) 



We can define an effective coupling strength, Q e ff, 
based on this attractor solution which would yield the 
same equation of state in the absence of the quartic term 



Qeff 



xmat{Q, Fq)- 



(36) 



The Einstein frame effective equation of state during the 
MAT-Q attractor is 



(37) 



C. Quartic term 

We now consider the implications of an additional 
quartic term f qU art<f> 4 m the action, with the coupling 
f quart parameterized as in (3). The Friedmann equation 
gives the fractional matter density, O m , in terms of these 
variables 



n m (a) 



I 



V 



iQF q x 



(31) 



While the non-minimal coupling, Q, in essence, speeds 
up the scalar's evolution, and increases we, the quartic 
term, F q , has the opposite effect, acting as a resistive 
force on the scalar and suppressing the kinetic fractional 
energy density in the MAT era , xmat- The effective 
equations of state in both the Einstein and Jordan frames 
are brought closer to w — as F q increases. In theory, 
therefore, this quartic coupling might lessen the tension 
between non-minimally coupled attractor solutions and 
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FIG. 1: The implications of the attractor solutions in the matter-dominated [left panel] and accelerative [right 
panel] eras when the leading order terms in the action are augmented by the quartic kinetic term parameterized by 
an amplitude F q . The panels show [top] the effective equations of state in the Jordan frame, u>j, and [middle] the 
Einstein frame, we along with [bottom] the effective coupling Q e ff [right] and potential exponent X e ff [left], 
relative to their fiducial values. Three values of non-minimal couplings ^i/2Q = 0.1, 0.3, 0.5 are considered. During 
the accelerative era Q only affects the conformal transformation and Jordan frame equation of state. 



constraints on the matter dominated era evolution from 
CMB distance measurements. 

In the accelerative era y 2 — 1 — x 2 — 36F q x 4 , and the 
attractor solution satisfies 

WqXACC + x \cc - ^ Xx acc = 0. (38) 

We can define an effective potential parameter, X e ff, that 
would give rise to the same dynamics in the absence of 
the quartic coupling, 

A e// = V6x ACC {F q ) (39) 
The Einstein frame effective equation of state 

1 



we = -l + ~\t ff (l + 4F q \i ff ) 



(40) 



As we see in Figure 1, the effect of an increasing quartic 
coupling is to reduce the effective equation of state at 
late times relative to the standard ACC-A value of w — 
— 1 + A 2 /3. We can see this consistently in the analytic 
solutions in the limit of small F q for the effective equation 
of state in the Einstein frame 



, A 2 4A 4 F, 
and in the Jordan frame 



3(1 - QX) V (1-QA) 



(41) 



(42) 



This should allow a larger range of values for A to be 
consistent with observations for F q > 0. 



D. Coupling to the Einstein tensor 

The presence of a direct coupling of the scalar to the 
Einstein tensor, with an amplitude parameterized by F c 
(4), modifies the Friedmann, acceleration and fluid equa- 
tions: 

fi m (a) = 1 - (1 + 18F c )x 2 - y 2 - z 2 , (43) 



i§{ 1+ T?k ll+im) * 2 ) = 1+ < 1+1SF c )x2 



and 



X = — ix — (1 + ^rr)X- 



1 



2 (1 + 6.F.) 



1 + 6F/ H 2 

[Xy 2 + QQ m (a)] . (45) 



Where \F C \ < 1/18 if the coupling to the Einstein tensor 
is to be subdominant to the canonical kinetic term. 
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While the new coupling does not introduce a new at- 
tractor, nor change the predictions for the RAD-A or the 
MAT- A attractor, it does modify the other two attractors 
from their nominal values, determined by the leading or- 
der terms. 

The MAT-Q solution now satisfies, 

= V6Q-3x-\/6Q:e 2 (1 + 24F c ) + (3 + 54F c )x 3 
+6\/6(l + 18F c )F c Qx 4 (46) 

leading to the solution 
1 



' = 1276^ ( - 3 + VITl6 ^ ) 



\Q (1 - 4F C Q 2 + ...) 



(47) 



the approximation holding for F c < 1, The effective 
equation of state in the Einstein frame is reduced for 
F c >0, 



we = ^ (-1 + V~l + 16F C Q2 



;Q 2 (1-4F c Q 2 + ...) 



as it also is in the Jordan frame, 

Wj = 



-l + 4Q 2 + ^l + 16FcQ 2 
3(l + 2Fc-2Q 2 ) 
4Q 2 / 4Q 2 F C 
3(1 - 2Q 2 ) ^ + 1 - 2Q 2 



(48) 

(49) 
(50) 



In the MAT— A solution, while the relative amplitudes 
of the kinetic and potential scalar densities are changed, 



1 



2(A-Q) 



y 



-Q(X-Q) (l + ^) 



A-Q 

the effective equation of state is unchanged 

Q 



W E 



A-Q' 



(51) 



(52) 



(53) 



The late time accelerative attractor equation, with 
y 2 = 1 - (1 + 18)a; 2 and z = 0, is 



-3x + 



V6 A 



3^6- 



2 l + QFc 1 + 6F C 

altering the ACC-A solution 

-1 - 6F C + y/l + 4F C (3 + 9Fc + A 2 ) 



Xx 2 (54) 



71 



1 - 6F r X 1 



2\/6F c A 
A 2 
6 



(55) 
(56) 



For small F c the effective accelerative equation of state 
is seen to become more negative for F c > 0, 



w E 



-1 + 



A 



2 r 



1 - 6F r 1 



A 2 
~6~ 



(57) 



and 



■l + M^(WJ^±?>| + ....(58) 



3(1 - QA) 



(1-QA) 



The numerical solutions are shown in Fig. 2. For increas- 
ingly positive values of F c the matter and accelerative 
era equations of state are lower, making them more con- 
sistent with data than F c = 0. While we show the impli- 
cations for < F c < 1 here, in our analysis to follow we 
will impose a restriction that \F C \ < 1/18 to ensure the 
term remains subdominant to those at leading order. 



E. A kinetic non-minimal coupling 

If we introduce a kinetic non-minimal coupling of the 
scalar field to matter, as in (5), the Friedman equation 
and acceleration equation in the Einstein frame arc un- 
changed, however the scalar equation of motion is modi- 
fied leading to a modified attractor equation 



H_ 

x +x H2 - 1 



x -3x + 



1 - 6F k x 2 



6F k x 2 
3 



(59) 



\ 2 I /^o O 'iFkin^'mX H 

2 V + V 2 ~ (l-6.F fc 2 2 )i^ 



This does not influence the effective equation of state in 
either the RAD-A or the RAD-null era, and the leading 
order matter and accelerative era attractors are unmod- 
ified by the inclusion of Moreover, the conformal 
transformation can still be calculated via (21) since only 
derivatives of fki n <P 2 = &F k x 2 enter the equation which 
are zero for attractor solutions. The coupling does, how- 
ever open up additional matter and accelerative era at- 
tractor solutions, with x = l/y/6F k . While we find that 
the numerical analyses do in some circumstances attempt 
to approach this attractor, the effect tends to be a tran- 
sitory, then returning to the leading order attractors. In 
its limit this attractor would lead to an ill-defined and 
unphysical = in the conformal transformation to the 
Jordan frame. Since the well-defined attractors are un- 
modified by this first order term, we do not investigate 
it further in the remainder of the paper. 



IV. COMPARISON WITH DATA 

A. Analysis Approach 

To investigate the impacts of the leading order and first 
order terms in the EFT, we numerically evolve the Ein- 
stein frame Friedmann, acceleration and scalar field fluid 




FIG. 2: The implications of the attractor solutions in the matter-dominated [left panel] and accelerative [right panel] 
eras when the leading order terms in the action are augmented by a coupling to the Einstein tensor parameterized 
by an amplitude F c . The panels show [top] the effective equations of state in the Jordan frame, wj, and [middle] the 

Einstein frame, we along with [bottom] the effective coupling Q e ff [right] and potential exponent X e ff [left], 
relative to their fiducial values. Three values of non-minimal couplings y/3/2Q =0.1, 0.3, 0.5 are considered. During 
the accelerative era Q only affects the conformal transformation and Jordan frame equation of state. 



equations, and simultaneously use the conformal trans- 
formation to calculate the Jordan frame variables. We 
compare the data to quantities in the Jordan frame since 
observations such as fluctuations in the CMB and the 
rcdshift of supernovae are reported assuming that the 
baryons are minimally coupled to gravity. We hence de- 
fine the present day epoch as aj = 1 and using the Jordan 
frame redshifts and distance measures as the physical ob- 
servables. 



package 3 . To ensure convergence we apply the Gelman 
and Rubin variance of chain mean/mean of chain vari- 
ances R statistic for each parameter on 8 or more chains. 
All MCMC runs have a convergence of R < 0.1 or lower. 

We consider constraints from the "Union 2.1" com- 
pilation 4 of 580 Type la supernovae observations [54], 
with redshifts in the range < z < 1.414. We compare 
predicted distance modulus estimates for the MCMC sce- 
narios for each supernovae, at redshift z, 



To establish the cosmological constraints on the EFT 
parameters, we perform a Monte Carlo Marcov Chain 
(MCMC) analysis assuming flat priors on Q m and 
H 0l along with leading order parameters Vo, A and Q. 
We consider constraints on the leading order terms plus 
each first order term separately to understand the indi- 
vidual effects of each. We assume flat priors on the quar- 
tic coupling, —10 < F q < 10, and the coupling to G^, 
< F c < 1/18 where the upper limit ensures sub domi- 
nance to the leading order terms. Scenarios with a Gauss 
Bonnet term are investigated with exponent < /i < 70 
and magnitude given by Fq — pgbFq 11 * where Fjf 6 * was 
defined in (30) and — 10 < log pgb < 0. 



n{z) = 5 log 



Dl(z) 
IMpc 



25, 



(60) 



where Dl is the luminosity distance, against the observa- 
tions using the compilation's covariance matrix including 
systematic errors. 

To investigate geometric constraints from the CMB, 
we use the WMAP-7 [8] results. The CMB is sensitive to 
two distance ratios to decoupling, through the position of 
the peaks and the acoustic oscillations [55]: the acoustic 
scale at decoupling, 



1 A (Z*) =(l + <)2^, 



(61) 



Our ID and 2D constraints are obtained after 
marginalizing over the remaining parameters using the 
programs included in the publicly available CosmoMC 



3 http://cosmologist.info 

4 http://supernova.lbl.gov/Union/ 
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and the 'shift parameter', 



n = 



H- 

11 n 



(l + z*)D A (z*), 



(62) 



where Da{z*) is the angular diameter distance, and 
r s (z st ) is the sound horizon, at the redshift to decoupling 

As discussed in [55] , the definition of 1Z suppresses the 
influence of radiation, dark energy or curvature on the 
Hubble parameter at decoupling H(z*) but is used by 
convention. We use the CMB data vector given in [8] 
WMAP 7 results 



(iA 
n 

\zj 



( 302.09\ 

1.725 
V1091.3/ 



with the inverse covariance matrix 
/2.305 29.698 



CoV CMB 



V 



-1.333 \ 
6825.270 -113.180 
3.414 j 



The redshift of decoupling, z„, is obtained , to percent 
accuracy, from the fitting formula [56] 



92 



.91 



0.560 
l + 21.1(O b /i 2 ) 1 - 81 ' 
0.0783(O^ 2 )- - 238 



(63) 



(64) 



1 + 39.5(n b h 2 ) 763 ' 

z* = 1048(1 + 0.00124(lM 2 )-°- 738 )(l + gi (n m h 2 ) 92 ). 

(65) 

where h = HQ/l00kms~ 1 Mpc~ 1 . In principle a non- 
minimal coupling, Q, will affect the the matter domi- 
nated expansion era and consequently the accuracy of 
the fitting function. We find that the redshift to decou- 
pling is only changed by 2 x 10~ 3 %, and the effects on 
I A and 1Z are smaller, so that the effect of the coupling on 
the accuracy fitting function in not a significant concern. 

Baryon Acoustic Oscillations (BAO) are the imprint 
of the sound horizon at last scattering as a characteristic 
scale in the clustering of matter. When observed at dif- 
ferent redshifts the characteristic scale can be used as a 
Standard Ruler to estimate cosmological distances. 

The acoustic scale along the line of sight encodes in- 
formation about the Hubble parameter H whereas the 
tangential component measures the angular diameter dis- 
tance Da- Current measurement accuracy is not suffi- 
cient to measure H and Da separately, so observational 
radial and tangential measurements are typically com- 
bined into an effective, averaged scale [57] defined as 



D v (z) 



D A (z) 2 (l + z) 2 cz 
H(z) 



1/3 



(66) 



The most accurate measurements of the BAO results to 
date come from 2dFGRS, SDSS DR7, WiggleZ and BOSS 
spectroscopic redshift surveys. The results presented in 
[14] are based on the spectroscopic SDSS DR7 sample, in- 
cluding both LRG and Main galaxy samples in combina- 
tion with the 2dFGRS survey. The ratio r s (zd r ag)/D v (z) 
is given for two redshifts, z — 0.2 and z — 0.35, 



r s {z drag )/D v {0.35) = 0.109715, 
r s (z drag )/D v (0.2) = 0.190533 

with the inverse covariance matrix 



(67) 
(68) 



Cov 



-l 

BAO-.SDSS 



/30124 -17227' 
y 86977 



Zdrag is the comoving sound horizon at the baryon drag 
epoch, baryon decouple from photons. If the ratio 
3/9{,/4p 7 = 1 at z, then the drag epoch and decoupling 
occur simultaneously. For typical, cosmological scenar- 
ios, however > Zd rag - We calculate the redshift Zd rag 
using the fitting formula from [58] 



Zd 

h 
b 2 



1291 



(n h 2 ) - 251 



[l + fe 1 (r> h /i 2 )"2],(69) 



T + 0.659(fi /io) - 828 
0.313(fi ^ 2 )°- 419 [l + 0.607(ft h 2 ) - 674 ], (70) 
0.238(ft /i 2 ) - 223 , (71) 



accurate to a few percent. 

As the distance scale D v is highly degenerate with 
Q m h 2 the WiggleZ survey introduced the acoustic pa- 
rameter A(z) [15], 



A{z) = D v {z)Jn m H 2 /cz. 



(72) 



The WiggleZ survey provides BAO measurements 
for three redshifts complementary to those from the 
SDSS/2dFGRS surveys: A(0.44) = 0.474, A(0.6) = 
0.442, A(0.73) = 0.424. The inverse covariance matrix 
is given by 



^ 0V BAO:WiggleZ 



/1040.3 -807.5 336.8 \ 

3720.3 -1551.9 
\ 2914.9 J 



Finally, we include the recent Baryon Oscillation Spec- 
troscopic Survey (BOSS) survey results, D v /r s — 13.67± 
0.22 at z = 0.57 [16]. 

In combination, the data sets have 589 — d degrees of 
freedom where d is the number of MCMC parameters. 
For KCDM d — 3, leading order d = 6 and all other 
model have d = 7. 



Model 


f2 m 


IQI 




ACDM 


r, 2 Q1 +0.031 
U.Z31_q q 28 








+Q (leading order) 


n OQ1+ 033 
u.zy±_ 030 


< 0.043 


< 1.15 


+F C 


n 2q2+ - 033 

U.Z»Z_ 029 


< 0.043 


< 1.40 


+F q 


n 2 Qo+0.033 
u.z»o_ 030 


< 0.044 


< 2.00 



TABLE II: Summary of the 95% confidence level 
constraints from the MCMC analysis for all scenarios 
except that including the Gauss-Bonnet term. The \ 2 

for all models are the same as for ACDM= 547.0 



B. Findings 

1. Constraints on the leading order, quartic kinetic and 
Einstein tensor terms 

In Table II we summarize the results of the MCMC 
analysis for the EFT models involving only the leading 
order terms, and those in which a quartic coupling or a 
coupling to the Einstein tensor are present, in comparison 
to ACDM. These models modify the attractor behaviors 
in both the matter and accelerative eras. The minimum 
X 2 is the same for each scenario, x 2 = 547.0, equivalent to 
that for ACDM, reflecting that in spite of including one 
or two extra parameters, these scenarios can recreate, but 
not improve upon, the predictions for ACDM. 

Scenarios with increasing magnitudes of coupling and 
exponential exponent come into tension with the data as 
they increase the predicted equation of state parameter 
in the matter era (for Q ^ 0) and the accelerative era 
(for A ^ 0) relative to the ACDM prediction. In the 
minimally coupled scenario, one would require A < V2 
to achieve any acceleration ( and A < \/2 to have w « 
— 1). The presence of a non-minimal coupling, Q < 1/A, 
creates a more negative equation of state in the Jordan 
frame, allowing a larger range of A, including potentially 
A > y/2, to be consistent with the data. 

In Figure 3 we show the combined 2D constraints aris- 
ing from CMB, BAO and SN for the coupling Q and the 
exponential exponent A. We find ID marginalized errors 
on the coupling of \Q\ < 0.026(0.043) and exponential 
potential, A < 0.80(1.15) at the 68% (95%) confidence 
level. In [39] a similar analysis was performed, but in 
the context of constraints arising from a non-minimal 
coupling purely to cold dark matter, for which the com- 
parison with observations is performed in the Einstein 
frame. They found similar constraints of |Q| < 0.055 
and exponential potential, A < 0.95 at the 95% confi- 
dence level. The similarity can be understood in terms 
of the small values of Q\ as the conformal transforma- 
tion tends towards unity, the Einstein and Jordan frame 
become comparable. 

As discussed in section III, while the quartic kinetic 
term and the coupling to the Einstein tensor don't lead 
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to new attractor solutions, they do alter the leading or- 
der attractor solutions during the matter dominated and 
accelerative eras. Their effects become less significant, 
however, as the coupling, Q, becomes smaller. We find 
that the constraints on Q are little changed by the in- 
clusion of these terms as its magnitude is already tightly 
constrained, to be very small, by the data. During the 
accelerative era these terms have a more pronounced ef- 
fect, and their enhancement of the accelerative equation 
of state allows a broader range of A to be consistent with 
the data, as given in Table II. 

While the effects of Q and A are sensitive to their mag- 
nitude, not sign, the quartic coupling term modification 
is sign dependent; negative values of F q and F c increase 
the effective equation of state during both the matter 
and accelerative eras. If the quartic term is included, we 
find the data provide a lower bound with F q > 0.22 at 
the 95% confidence level, as shown in Figure 4. For the 
coupling to the Einstein tensor, we find no significant dif- 
ference in best fit likelihoods in the range we investigated 
0<F C <1/18. 



2. Constraints on the Gauss-Bonnet term 

The Gauss-Bonnet term gives rise to a new accelera- 
tive attractor solution through creating a minimum in 
the effective scalar potential, produced by the exponen- 
tial potential and the GB term. This enables the expan- 
sion history to exit the scaling solution during the matter 
dominated era into a stable de-Sitter point. The GB term 
does not affect the RAD or MAT era attractors, so we 
consider its impact on a simple non-minimally coupled 
quintessence model, with Q = 0. 

In this scenario our present day epoch is in the transi- 
tion period between the matter and de-Sitter eras char- 
acterized by a trough [49] and plateau in the effective 
equation of state as shown in Figure 5. Prior to the 
transition the scalar follows the MAT-A attractor with 
w e ff — 0. At the transition epoch the positive gradient 
of GB term temporarily dominates the field's equation 
of motion. The trough in the effective equation of state 
is generated as the scalar slows and instantaneously be- 
comes static, then <j> changes sign, and the equation of 
state increases again. Following this the scalar proceeds 
on a slow evolution towards the static de-Sitter solution 
(in the distant future) . Both before and after the trough, 
w e ff is greater than in ACDM scenarios. 

We note that wfhile in the matter and de-Sitter eras the 
fractional energy density in the GB term, Qgb = V&xv, 
tends to zero, it is non-zero and can be significant in 
this brief transition between the two. During this epoch 
therefore we are close to, if not at the point, of no longer 
satisfying the sub dominance criteria used to formulate 
the EFT action. 

The depth of the trough is enhanced by a smaller value 
of A (giving a larger kinetic energy during the MAT-A 
phase), or steeper GB coupling, /j,. Enhancing the depth 
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FIG. 3: Joint 68% (dark shaded) and 95% (light shaded) constraints for the leading order action on the fractional 
matter density today, O m and the coupling of the scalar field to matter, Q, [left panel] and the scalar potential 

exponent, A [right panel]. 




FIG. 4: Joint 68% (dark shaded) and 95% (light shaded) constraints on the quartic kinetic coupling, F q and the 
coupling of the scalar field to matter, Q, [left panel] and the scalar potential exponent, A [right panel]. The quartic 
term has a significant effect on reducing the equation of state in the accelerative era, allowing a larger range of A to 

be consistent with the data. The evolution is sensitive to the sign of F q leading to a lower bound on the coupling. 



of trough also reduces the duration of the transient fea- 
ture, however, so that there is a sweet spot optimizing 
the trough's observational impact. The redshift position 
of the transient is degenerately sensitive to //, A and V 
through (30) and can be tuned through varying pgb- 

As summarized in Table III, we find that this distinc- 
tive equation of state profile is well constrained by the 
BAO and SN data, allowing us to place constraints on /i, 
A and pgb, and the contribution of the GB term in this 



effective action. 

In Figure 6 we show the tension between the con- 
straints from CMB+SN and CMB+BAO in terms of the 
exponents \i and A and fi m . We find the SN data have the 
strong preference for the GB trough to be located around 
z ~ 0.2 — 0.3, where previous principal component anal- 
yses of SN have shown the best measured modes in w e f / 
are peaked [59]. The accelerative era starts somewhat 
later than in ACDM and corresponds to a significantly 
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FIG. 5: The effective equation of state, w e ff, as a 
function of redshift, z,for the Gauss Bonnet model in 
comparison to a best fit ACDM scenario (black solid). 
[Top panel] Comparison of the best fit prediction for the 
Gauss-Bonnet model overall [blue, dashed] and for 
models with a low value of the exponent, /i < 18 (red 
dotted), showing how the depth of the "trough" in w e ff 
is dependent on fi. [Lower panel] Comparison of the 
overall best fit for CMB+SN (blue, dot- dashed), 
CMB+BAO (red, dashed dotted) and CMB+BAO+SN 

(black, dotted) showing the tension between the 
preferred evolution histories for the different datasets. 



FIG. 6: ID likelihood contours for the Gauss-Bonnet 

parameter [i and the potential exponent A for the 
different data sets: CMB+BAO+SN (black solid line), 
CMB+BAO (blue dashed line) and CMB+SN (red 
dotted line ). There is a clear tension between the 
supernova and the BAO data sets. As discussed in the 

text, through their effects on w e //, the data provide 
both upper and lower bounds on A and /i in this model. 



fit SN and BAO data. In this case, the Gauss-Bonnet 
model fits the data significantly worse than ACDM, with 
a difference of x 2 (GB)-x 2 (ACDM)=17. 



higher Q, m . Values of A ~ 4 — 5 are preferred to cre- 
ate a trough that is deep but still sufficiently broad to 
make the redshift averaged effective equation of state as 
consistent with that for ACDM as possible. By contrast 
the BAO place tight constraints on the matter density 
today. The average equation of state in the low red- 
shift regime is somewhat larger than that predicted in 
ACDM, which results in a comparatively higher Q m be- 
ing preferred to fit the BAO data. This is consistent with 
the standard degeneracies found in BAO constraints of a 
constant equation of state, e.g. [15]. The BAO are less 
sensitive to shape and depth of the trough than the SN, 
and place weaker constraints on A, /i and Pgb- 

A previous analysis of this scenario [53], only includ- 
ing the SDSS BAO survey, found reasonable agreement 
with the data when only considering CMB+BAO con- 
straints. We find that when we also include the WiggleZ 
and BOSS data the combined fit is not good and both 
the CMB+SN and CMB+BAO yield a worse fit than 
ACDM. The major issue arises when trying to jointly 



V. CONCLUSIONS 

We have considered the existence of dynamical attrac- 
tor solutions in a general approach to dark energy model 
building, utilizing methods of effective field theory. A 
wide range of dark energy and modified gravity models 
are able to be described using this approach, so that it 
forms a useful phenomenological link between underly- 
ing theories and observations. Dynamical attractors are 
powerful because of their weak sensitivity to initial con- 
ditions, which helps to alleviate some of the fine-tunings 
required to make models observationally viable. This ro- 
bustness can also be useful in determining constraints on 
the action of the effective theory and in the presence of 
stark tensions with astrophysical observations, can help 
isolate terms which are strongly disfavored - helping fo- 
cus model building efforts. 

We considered both the analytical and numerical pre- 
dictions for the cosmic expansion history and obtained 
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CMB+SN 
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u - za -0.03 


3.2 


n S5+ 03 






-2-9±i:3 


CMB+BAO+SN 
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TABLE III: Summary of the 95% confidence level constraints from Gauss-Bonnet MCMC analysis. The SN data 
provide the best constraints on GB parameters \i and A and pq b ■ There is significant tension between the 
constraints from SN and BAO datasets separately. The compromise, when the datasets are considered in 
combination, is a significantly worse fit to the data than ACDM. 
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FIG. 7: Joint 68% (dark shaded) and 95% (light 
shaded) constraints for the fraction matter density, VL m 

and for potential exponent A in the Gauss-Bonnet 
model for the different data sets: CMB+SN (red solid 
contours), CMB+BAO (yellow solid contours), and in 
combination CMB+BAO+SN (blue dashed lines). 
There CMB+SN and CMB+BAO data have almost 
orthogonal dependencies in this projection, highlighting 
the origins of the tension between them. 



numerical constraints on the effective theory in light of 
recent CMB, Type la SN and BAO constraints. In the 
Einstein frame attractor solutions exist that predict ob- 
servationally consistent radiation (RAD-null) and accel- 
erative eras (ACC-A), but require tight constraints on the 
presence of a non-minimal coupling to matter, Q, to give 
a viable matter era evolution (MAT-Q) . We have shown 
that the addition of terms which are quartic in the time 



derivative of the scalar field, scalar couplings to the Ein- 
stein tensor, and a Gauss-Bonnet term can all lead to 
modifications of the expansion history in the matter and 
accelerative eras. The quartic and Einstein tensor terms 
both modify the existing attractor solutions, creating ef- 
fective couplings and potential exponents Q e j / and A e / / 
that can reduce the effective equation of state parameter 
in both eras, and improve the fit to the data. They have 
a limited impact on the constraints of Q in the matter 
era, because their effect is diminished when Q is small. 
They can, however, have a significant affect on expand- 
ing the range of potentials consistent with the data by 
increasing the range of A allowed. 

The Gauss-Bonnet coupling opens up a new late-time 
de-Sitter solution (ACC-GB) induced by the creation of 
a minimum in the effective potential formed from V(<f>) 
and the GB term. In this scenario, our current epoch is 
a transitory era between matter domination and the de- 
Sitter phase, with a characteristic evolution in the equa- 
tion of state parameter that is constrained well by the 
data. In particular, a tension exists between the BAO 
and Type la supernovae, through their different redshifts 
sensitivity, when fitting the GB model to the data. In 
combination they rule out this scenario at large signif- 
icance, with a best fit x 2 {GB) - x 2 (ACDM)=17. This 
model, as well as being observationally disfavored, in- 
volves a transient epoch where the GB term becomes 
comparable with the leading order potential, so that on 
theoretical grounds one must also be cautious and ensure 
that the EFT remains valid. 

The complementarity of the SN, BAO and CMB dis- 
tance measures, has enabled us to place constraints on 
the cosmological background evolution by constraining 
the EFT couplings of higher dimensional operators. Our 
general approach is limited to perturbatively constructed 
backgrounds and does not apply to models where the 
background evolution becomes strongly coupled or highly 
non-linear. This means that we are unable to capture 
screening effects that could be an important additional 
signature for some models of modified gravity, as well 
as models such as k-essence where a perturbative de- 
scription of the background is not possible. To utilize 
the EFT approach in these models, one must put aside 
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the possibility of constraining the background evolution, 
and instead construct an effective theory for the per- 
turbations around an assumed ACDM-like background. 
Such an approach has already been taken for models 
of quintessence [31], where it was shown that instabil- 
ity issues - which are common place in such models - 
can be addressed within the context of the EFT. Work 
in preparation [60], will extend this approach in much 
the same spirit as here, but with an emphasis on con- 
straining the EFT of the perturbations. This will extend 
our analysis of cosmological constraints to the EFT of 
the perturbations, and their implications for complemen- 
tary constraints from large scale structure measurements 
of weak lensing, galaxy position and peculiar velocity 
fields. These correlations, comparing and contrasting rel- 
ativistic and non-relativistic tracers, could be a powerful 



probe of the broad range of gravitational modifications 
described by the EFT [61]. 
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